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Abstract 

Mass spectra and semileptonic decay rates of doubly heavy baryons 

■ are studied in the framework of the relativistic quark model in the 
quark-diquark approximation. 

p ■ The description of doubly heavy baryon properties acquires in the last 

years the status of actual physical problem which can be studied experi- 
mentally. The appearance of experimental data on B c mesons, heavy-light 
baryons stimulates the investigation of heavy quark bound systems and can 
^ ■ help in discriminating numerous quark models. Recently first experimen- 

tal indications of the existence of doubly charmed baryons were published by 
SELEX [1J. Although these data need further experimental confirmation and 
clarification it manifests that in the near future the masses and decay rates of 
doubly heavy baryons will be measured. This gives additional grounds for the 
theoretical investigation of the doubly heavy baryon properties. The success 
of the heavy quark effective theory (HQET) in predicting properties of the 
heavy- light qQ mesons (B and D) suggests to apply these methods to heavy- 
light baryons, too. The semileptonic decays of heavy hadrons present also 
an important tool for determining the elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. 

Doubly heavy baryons occupy a special position among existing baryons 
because they can be studied in the quark-diquark approximation and the 
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two-particle bound state methods can be applied. The two heavy quarks 
compose in this bound diquark system in the antitriplet colour state 

which serves as a localized colour source. The light quark q is orbiting around 
this heavy source at a distance much larger (~ than the source size 

(~ 2/rriQ), see Fig. ^ The estimates of the light quark velocity in these 
baryons show that its value is v/c ~ 0.7 — 0.8 and the light quark should be 
treated fully relativistically. Thus the doubly heavy baryons look effectively 
like a two-body bound system and strongly resemble the heavy-light B and 
D mesons. Then the HQET expansion in the inverse heavy diquark mass 
can be performed. The ground state baryons with two heavy quarks can be 
composed from a compact doubly heavy diquark of spin or 
quark. According to the Pauli principle the diquarks (bb) or 
spin 1 whereas diquark (be) can have both the spin and 1. 

Here we study mass spectra and semi- 
leptonic decay rates of doubly heavy baryons 
using the relativistic quark model in the quark 
-diquark approximation. 

In the quark-diquark picture of doubly 
heavy baryons the bound states of two heavy 
quarks and of the light quark and the heavy 
diquark are described by the diquark wave 
function (^d) and by the baryon wave func- 
tion (^b), respectively. These wave functions 
satisfy the two-particle quasipotential equa- 
tion of the Schrodinger type [2] 




Figure 1: Schematic picture 
of doubly heavy baryon. 
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where the relativistic reduced mass is 
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and the center of mass energies of particles on the mass shell E\, E 2 are given 
by 

M 2 — ml + mj ^ M 2 — m 2 + m 2 
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Here M = Ei + E 2 is the bound state mass (diquark or baryon), m^ 2 are the 
masses of heavy quarks (Qi and Q 2 ) which form the diquark or of the heavy 
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diquark (d) and light quark (q) which form the doubly heavy baryon (B), 
and p is their relative momentum. In the center of mass system the relative 
momentum squared on mass shell reads 

2 [M 2 -(m 1 + m 2 ) 2 ][M 2 -(m 1 -m 2 ) 2 ] 

4M 2 • KV 

The kernel V(p, q; M) in Eq. (JTJ) is the quasipotential operator of the 
quark-quark or quark-diquark interaction. It is constructed with the help 
of the off-mass-shell scattering amplitude, projected onto the positive energy 
states. Here we closely follow the similar construction of the quark- ant iquark 
interaction in heavy mesons which were extensively studied in our relativistic 
quark model H] . For the quark-quark interaction in a diquark we use the 
relation Vqq = Vqq/2 arising under the assumption about the octet structure 
of the interaction. The quasipotential of the quark- ant iquark interaction Vqq 
is the sum of the usual one-gluon exchange term and the confining part which 
is the mixture of long-range vector and scalar linear potentials, where the 
vector confining potential contains the Pauli term. The explicit expressions 
for Vqq and Vj g and the details of the mass spectrum calculation are given 
in Ref. [5]. The quark masses have the following values m& = 4.88 GeV, 
m c = 1.55 GeV, m s = 0.50 GeV, m Ujd = 0.33 GeV. 

The masses of the ground state axial vector diquarks were found to be 
M£ v = 3.226 GeV, M^ v = 9.778 GeV, M£ v = 6.526 GeV, and the mass 
of the scalar diquark M^ c = 6.519 GeV. We calculated the mass spectra of 
light-quark - heavy-diquark system first in the infinitely heavy diquark mass 
limit and then with the inclusion of 1/Md corrections. In the infinitely heavy 
diquark mass limit its mass and spin decouple and the dynamics of the heavy 
hadron is determined by the light quark alone in accord with the heavy quark 
symmetry. Thus the properties of heavy-light mesons and doubly heavy 
baryons are similar in the limit tuq — > oo. Inclusion of the first order I/tuq 
corrections breaks the heavy quark symmetry and leads to different splittings 
in mesons and baryons jHlE]- Note that in our calculations we treated the 
light quark completely relativistically without applying unjustified expansion 
in inverse powers of its mass. 

The mass level orderings of S cc and E&& baryons are schematically shown 
in Fig. |21 There we first show our predictions for spectra in the limit when 
all 1/Md corrections are neglected. In this limit the P-wave excitations of 
the light quark are inverted. This means that the mass of the state with 
higher light quark angular momentum j = 3/2 is smaller than the mass of 
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Figure 2: Masses of S cc (left) and baryons (right) (in GeV). 

the state with lower angular momentum j = 1/2. Next we switch on 1/Md 
corrections. This results in splitting of the degenerate states and mixing of 
states with different j, which have the same total angular momentum J and 
parity. The fine splitting of P-levels turns out to be of the same order of 
magnitude as the gap between j = 1/2 and j = 3/2 degenerate multiplets in 
the Md — > oo limit. The inclusion of 1/Md corrections leads also to relative 
shifts of the hadron levels further decreasing this gap. As a result, some of 
the P-levels from different (initially degenerate) multiplets overlap; however, 
the centers of levels averaged over the heavy diquark spin remain inverted. 

s w 





Figure 3: Weak transition matrix el- Figure 4: The leading order contribu- 
ement of the doubly heavy baryon in tion to the diquark vertex func- 
the quark-diquark approximation. tion T. 



The consideration of semileptonic decays of doubly heavy baryons (bbq) 
or (bcq) to doubly heavy baryons (bcq) or (ccq) can be divided into two steps 
(see Fig. 3). The first step is the study of form factors of the weak transition 
between initial and final doubly heavy diquarks. The second one consists in 
the inclusion of the light quark in order to compose a baryon with spin 1/2 
or 3/2. 

In the relativistic quark model the transition matrix element between two 
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diquark states is determined jB] by the convolution of the wave functions ^ 
of the initial and final diquarks with the two particle vertex function T 



(d>(Q)\jW\d(P)) 



d 3 p d 3 q 
(2nY 



* < ,/ )0 (q)r M (p,q)* di p(p). 



(5) 



Here P = vMi, Mi, v denote the four- momentum, mass and four- velocity of 
the initial diquark Q b Q s and Q = v'Mf, Mf, v' denote the four- momentum, 
mass and four- velocity of the final diquark Q a Q s ; p and q are the relative 
quark-quark momenta. 

The leading contribution to the vertex function comes from the di- 
agram in Fig. 4 (index b denotes the initial active quark, index a the final 
active quark and index s the spectator) 

r M (p, q) = rW = u a (q x )7 M (l - 75K(PiK(q2K(P2)(27r) 3 5(p 2 - q 2 ), (6) 

where u(p) is the Dirac spinor. 

The transformation of the bound state wave function from the rest frame 
to the moving one with four-momentum P is given by |Hj 



(7) 



where R w is the Wigner rotation, Lp is the Lorentz boost from the diquark 
rest frame to a moving one, and D 1 ^ 2 (R) is the rotation matrix. 

Using this relation and the properties of the Dirac spinors and rotation 
matrices we can express the matrix element (jSj) in the form of the trace over 
spinor indices of both particles [7j. The final covariant expression for the 
transition matrix element reads 



(d'(Q)\J?\d(P)) fd 3 pd\ 

(27T) 
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l -Tr{^ df (Q,q) lll (l- l5 )^ d (P,p)}6 3 (p 2 -q 2 ), 
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where the amplitudes for the scalar (5") and axial vector (AV) diquarks (d) 
are given by 
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Here <&d(p) = ^d,o(p)/V2^d is the diquark wave function in the rest frame 
normalized to unity and the four- vector p = Lp(0, p). The argument of the 
5-function in Eq. (jHJ can be rewritten as 



P2 - q2 = q - p 



ts(p) + e s (q) (V _ y) 
w + 1 y h 



where w = (f • v'). The spectator quark contribution factors out in all decay 
matrix elements. They have a common factor 
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The weak current matrix elements, e.g., for the scalar to axial vector 
diquark transition (be — > cc) have the following covariant decomposition 



(AV(v',e')\jv\S(v)) ^ 

- = ih v {w)e tia ^e v p v^, 
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These transition form factors are expressed through the overlap integrals of 
the diquark wave functions and are given in Ref . [7] . These exact expressions 
for diquark form factors were obtained without any assumptions about the 
spectator and active quark masses. For the heavy diquark system we can ap- 
ply the v/c expansion. Then in the nonrelativistic limit we get the following 
expressions for the form factors 

h v (w) = [l + (w + l)f(w)]F(w), 
h Al (w) = h Aa (w) = [l + (w-l)f(w)]F(w), 
h A2 (w) = -2f(w)F(w), (15) 



where 



F(w) 



1/2 



W[W+1) \ J m 2 + ( w 2 _ D m 2 



and 



x /(S5^( p + Jri (v '- v) >*< (p) (16) 



/M = = • (17) 



The appearance of the terms proportional to the function f(w) is the result 
of the account of the spectator quark recoil. Their contribution is important 
and distinguishes our approach from the previous considerations [HIE]- We 
plot the function F(w) for bb — > be and be — > cc diquark transitions in Fig. E| 
The function F(w) falls off rather rapidly, especially for the bb — > be diquark 
transition where the spectator quark is the b quark. Such a decrease is the 
consequence of the large mass of the spectator quark and the high recoil 
momentum (g max ~ m& — m c ~ 3.33 GeV) transfered. 

The second step in studying weak transitions of doubly heavy baryons 
is the inclusion of the spectator light quark in the consideration. We carry 
out all further calculations in the limit of an infinitely heavy diquark, — > 
oo, treating the light quark relativistically. The transition matrix element 
between doubly heavy baryon states in the quark-diquark approximation (see 
Figs. 3 and 4) is given by [cf. Eqs. © and (jSJ)] 



B'{Q)\JW\B{P)) rd^pSq 



^,Q(q)(rf / (Q)|Jr|rf(P))^B,p(p)5 3 (p g -q (i 



(18) 
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Figure 5: The function F(w) for the bb — > 6c (left) and 6c — ► cc (right) quark 
transitions. 

where ^b,p(p) is the doubly heavy baryon wave function; p and q are the 
relative quark-diquark momenta. The baryon ground-state wave function 
*&B,p(p) is a product of the spin- independent part ^b(p) satisfying the related 
quasipotential equation (JTJ) and the spin part Ub(v) 



The baryon spin wave function Ub is constructed from the Dirac spinor 
u q {y) of the light spectator quark and the diquark spin wave function. The 
ground state spin 1/2 baryons can contain either the scalar or axial vector 
diquark. The former baryon is denoted by Sqq, and the latter one by Hqq/. 
The ground state spin 3/2 baryon can be formed only from the axial vector 
diquark and is denoted by Hqq/. 

The amplitude for the Hqq s — > ^q'q s transition in the infinitely heavy 
diquark limit is given by the following expression 



where rj(w) is the heavy diquark - light quark Isgur-Wise function which is 
determined by the dynamics of the light spectator quark q 



*b,p(p) = V B (p)U B (v). 



(19) 





(21) 
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We plot the Isgur-Wise function r](w) in Fig. |H1 In the nonrelativistic limit 
for heavy quarks the diquark form factors hi(w) contain the common factor 
F(w)rj(w). 

Our results for the semileptonic de- x r — 

cay rates of doubly heavy baryons E bb " : 

and Sft c are compared with previous pre- 
dictions in Tabled The results of differ- °' 6 
ent approaches differ substantially. Most 0.4 
of previous papers (TU1 HH UJ] give their 2 
predictions only for selected decay rates. 

Their values agree with our in the or- 1 1 - 05 1,1 1,15 1,2 
der of magnitude. Our predictions are pigure g. Thg Isgur . Wise func . 
smaller than the QCD sum rule results ^ ^ q{ ^ ^ quark _ 
dH by a factor of ~ 2. This can be a re- heayy diqu&rk bound sygtem 
suit of our treatment of the heavy spec- 
tator quark recoil in the heavy diquark. On the other hand the authors of 
Ref. [12] using for calculations the Bethe-Salpeter equation give more decay 
channels. Their results are substantially higher than ours, for some decays 
the difference reaches almost two orders of magnitude which seems quite 
strange. E.g., for the sum of the semileptonic decays E bb — > 'E b '* Ref. [T2] 
predicts ~ 6 x 1CT 13 GeV which almost saturates the estimate of the total 
decay rate r|°* al ~ (8.3 ± 0.3) x KT 13 GeV [TH] and thus is unlikely. 
This work was supported in part by the DFG grant No. Eb 139/2-2. 
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